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We present a search for Kaluza-Klein (KK) particles predicted by models with universal extra 
dimensions (UED) using a data set corresponding to an integrated luminosity of 7.3 fb~^, collected 
by the DO detector at a pp center of mass energy of 1.96 TeV. The decay chain of KK particles can 
lead to a final state with two muons of the same charge. This signature is used to set a lower limit 
on the compactification scale of R^^ > 260 GeV in a minimal UED model. 

PACS numbers: 04.50.-h, 14.80.Rt, 13.85.Rm 



The existence of extra dimensions in addition to the 
3 + 1 dimensions of space-time has been postulated as 
a possible solution to the problem of the large hierarchy 
of scales in the standard model (SM). In models with 
universal extra dimensions (UED) all particles propagate 
in the extra dimensions In this Letter, we study a 
minimal UED (mUED) model, which has only one extra 
dimension [2]. 

Each SM particle in the mUED model is associated 
with a set of excited Kaluza-Klein (KK) states when 
viewed in 3 -f 1 dimensions. Since compactification of the 
extra dimensions leads to periodic boundary conditions, 
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the KK states have discrete masses of the order of R~^, 
where R is the radius of the compact dimension. If one- 
loop corrections are applied, the mass spectrum of the 
KK modes also depends on a cutoff scale for boundary 
terms, which is chosen to be 10 TeV [2] . Gluon KK modes 
[gi) are the heaviest particles, followed by quarks (SU(2) 
doublet Qi or singlet qi), gauge bosons [Zi/Wi), lep- 
tons (SU(2) doublet Li or singlet £i) and the KK photon 
(71), which is the lightest KK particle (LKP) and does 
not decay. The LKP is also a dark matter candidate j^]. 

Previous searches for KK particles predicted by a mod- 
ified UED model have been conducted by the DO Col- 
laboration at the Fermilab Tevatron Collider and by the 
ATLAS [H^l Collaboration at the CERN Large Hadron 
Collider (LHC) in diphoton final states. In this model, 
gravity mediation allows decays of the LKP into a pho- 
ton plus a light KK graviton. A search for events with 
two leptons of the same charge, which can be interpreted 
in a mUED model, has been performed by ATLAS @, 
but no dedicated study of the mUED model has so far 
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been performed at a collider. Constraints from preci- 
sion electroweak data [ij , from the measurements of the 
muon anomalous magnetic moment and of the 6 — ^ 57 
branching ratio Q, indicate that the scale can be 
as low as ~ 300 GeV, making KK particles accessible at 
the Tevatron. In this Letter, the first search for mUED 
in the final state with two leptons of the same charge 
is presented using data corresponding to 7.3 fb~^ of in- 
tegrated luminosity collected by the DO detector [§] at 
Vs = 1.96 TeV. 

At the Tevatron, KK gluons or quarks are mainly pro- 
duced in pairs, as shown in Fig. [TJ In the subsequent 
cascade decay, up to four charged leptons are produced. 
Since the masses of the extra particles predicted by the 
mUED models are nearly degenerate, the leptons are 
emitted with low transverse momentum and might es- 
cape detection. In this analysis, we select events with 
two muons of the same charge. 




FIG. 1: (a) Production of a pair of KK quarks (QiQi). (b) 
Decay of a KK quark into a jet, two oppositely charged lep- 
tons, and the LKP. Double lines indicate KK excitations. A 
similar cascade decay occurs for the second KK quark leading 
to several leptons of the same charge in the final state. 

The DO detector [§] consists of tracking systems and 
calorimeters. The innermost part is a tracking system 
where charged particles are detected by the silicon mi- 
crostrip (SMT) and central fiber (CFT) tracking detec- 
tors, located within a 2 T solenoid. The tracking sys- 
tem is surrounded by a liquid- argon/uranium calorime- 
ter. Particle energies are measured in the electromagnetic 
and hadronic calorimeters within a pseudorapidity range 
of \r]\ < 4.2 0. Jets are reconstructed with a cone algo- 
rithm using a radius of 7?. = 0.5 [lo| in the calorimeter. 
The central and forward muon detectors are composed 
of a layer of wire chambers and scintillators in front of 
a 1.8 T toroid magnet and two layers outside the toroid. 
Missing transverse energy is measured from the vector 
sum of the calorimeter cell energies in the xy plane 11 1 . 
A correction for the energy response of muons, electrons 
and jets is applied. 

The backgrounds from Z-|-jets, VF-Mets, and ti pro- 
duction are modeled by the ALPGEN 
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Monte Carlo 

(MC) event generator, interfaced with PYTHIA [l^ for 
showering and hadronization. Diboson production ( WW, 
WZ and ZZ) is simulated by PYTHIA. The CTEQ6L1 
parametrization of the parton distribution functions 
(PDF) are used [l3]- Higher order cross sections for 
diboson and W/Z +}ets production are calculated by 
MCFM [1^ and the cross section of tt pair production 



is taken from 16[ . The dominant source of background 
is pairs of muons with the same charge from heavy fla- 
vor jets. This background contribution is estimated from 
data and is described in detail in the following. 

Signal MC events are generated for 9 different val- 
ues of covering the range from 200 to 320 GeV in 
steps of 15 GeV, using pythia with the CTEQ5L PDF 
parametrization p/7j. The production cross sections and 
masses of KK particles are taken from PYTHIA. They are 
given in Table H] for each and with all KK gluon and 
quark production modes included. All decay mechanisms 
leading to like-charge dimuon final states are taken into 
account. The decay branching fractions for all KK par- 
ticles are given by the mUED model. After simulating 
all cascade decays, approximately 1% of events have two 
like-charged muons. 



TABLE I: Masses of KK particles for each R ^ value used in 
the MC generation with corresponding total production cross 
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Masses (GeV) 




Cross Section 


(GeV) 


71 


Zi 


91 
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Qi 


(pb) 


200 


201 


230 


269 


207 


249 


34.9 ± 0.2 


215 


216 


245 


287 


222 


266 


20.4 ±0.1 


230 


231 


260 


305 


238 


283 


12.1 ± 0.1 


245 


246 


274 


323 


253 


300 


7.24 ± 0.05 


260 


261 


289 


341 


268 


317 


4.39 ± 0.03 


275 


276 


304 


359 


284 


334 


2.69 ±0.02 


290 


291 


319 


377 


299 


351 


1.65 ±0.01 


305 


306 


335 


395 


314 


368 


1.02 ±0.06 


320 


321 


350 


413 


330 


385 


0.63 ±0.01 



All signal and background MC events pass through the 
full GEANT-based simulation of the detector [l^l and are 
reconstructed using the same algorithms as used for data. 
To simulate detector noise and multiple pp interactions, 
MC events are overlaid with data events from random 
beam crossings. 

Events are selected by requiring that they pass at least 
one single muon trigger condition. We require that each 
event must contain at least two muons of the same charge. 
The track in the muon system must be matched to a track 
in the central tracking system with detector 1 77 1 < 1.5. We 
reject cosmic rays by requiring the associated scintillator 
hits in the muon system to be consistent with originating 
from a pp collision, the distance of closest approach (dca) 
of the muon tracks to the pp interaction vertex to be less 
than 0.05 cm, and the differences between the z coor- 
dinates of the dca of each muon and the pp interaction 
vertex to be < 1 cm. We require A(/)^^ < 2.9 rad for the 
azimuthal angle between the two muons in order to reject 
multijet background in which there are back-to-back jets 
which can contain muons from semileptonic meson decay, 
and A0^^ > 0.25 rad to reject muons originating from 
the same jet. The transverse momenta of the leading and 
next-to-leading muons have to be 15 < pTi < 200 GeV 
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and > 10 GeV, respectively. The invariant mass of 
the muon pair must be in the range M^^ < 250 GeV. 
The upper hmits on the muon px and dimuon invariant 
mass reduce the number of events with one of the muon 
charges mis-reconstructed. In addition, we reject events 
with < 25 GeV, since multijet background dominates 
at low ■ 

To discriminate between isolated muons from signal 
and muons contained in jets, we define the isolation in 
the calorimeter, 1'="', as the sum of the energy deposited 
in the calorimeter within an annulus of 0.1 < TZ < 0.4, 
divided by the muon px, and the isolation in the track- 
ing detector, I*'"'^, as the sum of the pt of all charged 
particles within a cone oi TZ = 0.5 around the muon, ex- 
cluding the muon itself, divided by pt^. At least one of 
the two muons is required to be isolated with Z"^"^' < 0.4 
and < 0.12. 

We estimate the multijet background by defining a sig- 
nal and a background enriched sample. The signal en- 
riched sample comprises events where the requirement on 
the isolation of the second muon in the tracking detec- 
tor is relaxed to X*'''^ < 0.25. To define a background 
enriched sample, we require that the second muon fails 
the isolation requirement. A normalization factor is cal- 
culated for each jet multiplicity, given by the ratio of the 
number of events in the signal and background enriched 
samples in a multijet dominated region where the pt of 
the most isolated muon is 5 < pt < 10 GeV. The multijet 
background is determined by multiplying the background 
enriched sample with this normalisation factor in the re- 
gion where the pr of the most isolated muon is > 10 GeV. 

In the high px region, there is a significant contribution 
to the background enriched sample from SM processes 
other than multijet, particularly from I1^-|-jets produc- 
tion, with an isolated muon from the W boson decay 
and a non-isolated muon embedded in the jet. This SM 
background is modeled by applying the selection used to 
define the background enriched sample to MC events. 
The resulting SM background distributions, corrected 
with the same normalization factors as used for data, are 
subtracted from the corrected distribution of the back- 
ground enriched data sample to obtain an estimate of the 
multijet background. 

Events from Z — >■ fj,fj, decays enter our sample mainly 
through mis-reconstruction of one of the muon charges. 
This background is estimated from simulation. We use 
data to determine the corresponding systematic uncer- 
tainty by comparing the two independent charge mea- 
surements in the central tracking detector and the muon 
detector [l^. We determine the rate of muons with an 
incorrectly measured charge by counting the number of 
events where the two measurements disagree. The num- 
ber of charge-flip events in data is in good agreement with 
the number of like-charged Z — > /i/i events predicted by 
the simulation. 

To improve the discrimination between signal and 
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(b) D0, 7.3 fb"' 
FT^ = 260 GeV 
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FIG. 2: Distribution of (a) and (b) pti for data and 
background, compared to a signal with = 260 GeV. The 
shaded band shows the statistical uncertainty on the back- 
ground estimation. All entries exceeding the range of the 
histogram are added to the last bin. 



background a boosted decision tree (BDT) [20] is 
trained using pT^^, pjb, -^/i/^, £ind the num- 

ber of jets as input variables. Further input vari- 
ables are the of the fit of the muon tracks to re- 
ject badly reconstructed muons, the transverse masses 



Mt — y2_^tPt(1 — cosA0(^t,m)) calculated sepa- 
rately for each muon, and the scalar product of and 
PT2- This rejects multijet background events with small 
or • 

The distributions of and pxi for data and back- 
ground are shown Fig. [51 together with signal for = 
260 GeV. The BDT is trained separately for each signal 
point with different R~^. The BDT output distribution 
is shown in Fig. [31 and the event yields in data are com- 
pared to the expected number of background and signal 
events in Table [ill Signal is concentrated in the region 
where the BDT output is > 0. 

Systematic uncertainties on the normalisation of both 
background and signal, including their correlations, are 
taken into account. These include theoretical uncer- 
tainties on SM background cross sections (7% — 15%) 
and the uncertainties on the rate of multijet background 
(40%), signal efficiency from choice of PDF parametriza- 
tion (4%), integrated luminosity (6.1%), jet energy scale 
(4%), rate of Z/"f* production with an incorrectly recon- 
structed muon charge (21%), trigger efficiency (6%) and 
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TABLE II: Expected number of events for backgrounds, event 
yields in data and expected number of events for a signal with 
— 260 GeV after the final selection and requiring a BDT 
output > (for illustrative purpose, this cut is not used for 
limit estimation). The total uncertainties are also given. 



Process 


Final Selection BDT output > 


Diboson 


21 ±3 


6± 1 


Z+jets 


39 ±9 


13 ±3 


W+jets 


109 ± 14 


38 ±5 


tt 


6± 1 


2± 1 


Multijet 


95 ±41 


63 ±27 


Total Background 


271 ± 45 


123 ± 28 


Data 


273 


126 


Signal 


18 ± 1 


18 ± 1 



« 10^ 
c 

> 

UJ 
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-t- Data 
I I Multijet 
□ Z + jets 

■ W + jets 

■ Diboson 

— Signal 



fi"^ = 260 GeV"^ D0, 7.3 fb"^ 




0.5 
BDT output 



i 0.25 
m 0.2 

X 

0.15 
0.1 
0.05 










D0 Run II, L = 7.3 fb ' 




— ^ Observed Limit 




Expected Limit 




± 1 SD Expected 




1 1 ± 2 SD Expected 




Theoretical 





200 



220 



240 



260 



280 300 
R"' [GeV] 



FIG. 4: Observed and median expected 95% CL limits on 
aB{fi'^ fi^) as a function of compared to aB{^^ cal- 
culated with the mUED model. The bands represent ±1 and 
±2 standard deviations around the median expected limits. 



dia); Colciencias (Colombia); CONACyT (Mexico); NRF 
(Korea); CONICET and UBACyT (Argentina); FOM 
(The Netherlands); STFC and the Royal Society (United 
Kingdom); MSMT and GACR (Czech Republic); BMBF 
and DFG (Germany); SFI (Ireland); The Swedish Re- 
search Council (Sweden); and CAS and CNSF (China). 



FIG. 3: Distribution of the BDT output for data and back- 
ground, compared to a sig nal with R'^ = 260 GeV. The 
shaded band shows the statistical uncertainty on the back- 
ground estimation. 



muon reconstruction and isolation (2%). 

The entire differential distribution of the BDT output 
is used to set limits on the product of the cross section 
and the branching ratio with the CLg method 2l| and 
a profiling technique to reduce the impact of systematic 
uncertainties [l^l . Observed and median expected upper 
limits at 95% CL on the product of the cross section 
a and the branching fraction B into like-charged muon 
pairs are shown in Fig. U as a function of . The 
theoretical cross section of the mUED model intersects 
the expected limit at R^^ — 275 GeV and the observed at 
R^^ = 260 GeV, which corresponds to a mass of 317 GeV 
of the lightest KK quark in the mUED model. 

In summary, we present a search for extra dimensions 
in the mUED model using 7.3 fb^^ of integrated luminos- 
ity collected by the DO experiment. The first direct lower 
limit on the compactification scale of the extra dimension 
in the mUED model is set as R'^ = 260 GeV. 
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